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Abstract — Several robust eigenstructure assignment meth- 
ods for MIMO systems try to exploit non-uniqueness of 
controller matrix to optimize certain performance index. In 
general these methods lead to computational intractability 
or become too conservative design. Therefore the scope of 
using randomized algorithms as an alternative exists. Based 
on well established statistical learning theory algorithms this 
paper provides methodologies to design robust controller under 
structured uncertainty by eigenstructure assignment method. 
Also a Monte Carlo based computational procedure for a low 
sensitivity regional pole assignment problem is illustrated. With 
randomized methods, a controller can be easily found working 
satisfactorily for most of the time, while the price to pay for is 
its probabilistic completeness. 

I. INTRODUCTION 

Eigenstructure assignment is a well established method 
that directly incorporates classical time domain design spec- 
ifications. Initiated by Wonham [1] and later, Moore [2] 
described the non-uniqueness of the controller which can 
assign a self conjugate set of eigenvalues for MIMO system. 
The seminal paper of Srinathkumar [3] clearly describes 
the freedom available to the designer in selecting maxi- 
mum number of eigenvalues and number of entries in each 
corresponding eigenvector. There exist several iterative and 
non-iterative methods which exploits this freedom to design 
appropriate controller for meeting the design specifications. 
Like, since the early design periods, attempts have been made 
to assign eigenvalues having low sensitivity with respect to 
perturbation in system matrices for better robustness. These 
can be thought of some what indirect approaches to design 
a robust controller. In these methods, the sensitivity measure 
considered in general, is the overall eigenvalue condition 
number which is a function of quadratic norm of modal 
matrix. There exists several methods for the problem of 
exact eigenvalue placement with minimum overall eigenvalue 
condition number. One can find in [4] details of low sensitiv- 
ity eigenstructure assignment methods for exact eigenvalue 
location. In practice, it is hardly required to place the 
eigenvalue exactly. Instead, often the design requirement is to 
place the poles in pre-specified regions in the complex plane. 
A low sensitivity eigenvalue assignment within a specific 
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region is considered by Lam and Tam [5]. Their approach is 
based on gradient based minimization technique to find the 
optimum pole location, which corresponds to a minimum 
Frobenius condition number of the closed loop state matrix. 
In general such a design is over conservative when com- 
pared to designs based on minimization of quadratic overall 
eigenvalue condition number. Note that quadratic condition 
number minimization with respect to a relaxed pole location 
is in general a non convex problem. 

In contrast to minimizing overall eigenvalue sensitivity, 
several direct attempts have been made to ensure both 
time domain and frequency domain robust stability and 
performance measures. This is done by combining frequency 
domain stability and performance measures along with time 
domain design specification to realize joint optimum robust- 
ness. In [4] and [6], minimization of an infinity norm of 
mixed sensitivity measure over all possible controller is done 
using genetic algorithm and gradient based methods. But 
as shown by researchers, //«, norm minimization methods 
results in very conservative controller design. 

From early 1990's control researchers began to consider 
the computational complexity of several control problems. 
One can find a survey of existing results regarding computa- 
tional complexity of many robust control problems in [7]. In 
particular, simultaneous stabilization of linear systems with 
static state/output feedback when the controller parameters 
are constrained is proved to be NP-Hard. (Clearly one can 
relate the problem of simultaneous stabilization of linear 
systems by constrained controller to the problem of robust 
constrained eigenstructure assignment methods, since in the 
later case, the designer generally put constraints on some 
entries of desired eigenvector matrix to shape the time 
response/ mode decoupling). Details can be found in [8] and 
[9]. 

With such pessimistic result, an alternative approach that 
is recently gaining popularity is to use randomized algo- 
rithm to break the curse of dimensionality. In general these 
algorithm are not required to work in all cases, but rather 
most of the cases. These algorithms have polynomial time 
complexity. The price one pays is that, these algorithms 
are only probabilistically complete. Early attempt to bring 
randomization into control is initiated by Ray and Stengel 
[10]. More recently Khargonekar and Tikku [11], Tempo and 
Dabbane [12] and Vidyasagar [13] contributed to this field. 

With this motivation, in this paper a randomized approach 
to low sensitivity eigenstructure assignment for relaxed 
pole location is developed. Also, application of random- 
ized algorithm originated from statistical learning theory 



on eigenstructure assignment in designing robust control 
under structured/parametric uncertainty is shown. The paper 
is organized as follows: in Section II classical eigenstructure 
assignment method is reviewed. In Section III a robustness 
measure is defined and consequently general nature the 
objective functions is investigated. A randomized algorithm 
method to solve low sensitivity eigenstructure assignment 
problem is presented. Section IV contains a direct application 
of statistical learning theory to the design of constrained 
robust controllers against parametric uncertainty, while in 
subsequent section we discuss the result and conclude. 

II. Eigenstructure Assignment Method for LTI 

Systems 
Consider a MIMO linear time invariant system described 
by 

x = Ax + Bu (1) 

y = Cx (2) 

where the state vector i£l", control input vector u G R m , 
output vector y e W and A E R" x ", B e K" xm and C e W x " 
system matrices. Assume rank(B) = m and rank(C) = r 
(This full rank assumption is not strictly required). Define 
a required closed-loop self conjugate eigenvalue set A = 
{A; : A,- € C, i = 1,2, ...«}. For simplicity assume all required 
eigenvalues are distinct. The linear output feedback law is 
applied to the open loop system as u = Ky where K G K mx '' 
such that the closed loop eigenvalue specified in A is attained. 
Let us denote the i ,h left and right eigenvector corresponding 
to Xj as Lj and /?, respectively. The modal and spectral 
matrices defined as 



R = [Ri R 2 ...R„]eR" xn 
L = [L, L 2 ...L„]eR nxn 

and D A = diag[X u Xz,~X] £ K" x ". 

Then the forced response of the system is given by 



(3) 
(4) 



III. Low Sensitivity Eigenstructure Assignment 
under Regional Pole Constraints 

Define A = {Ai,A2,...A„} to be a self conjugate set of 
eigenvalues. Now given a system matrix (1) one can always 
select a K = WR^ 1 such that A is attainable. Define J(R) = 
rj(R) = \\R\\ 2 \\R~ l \\ 2 as the quadratic condition number. 
The control problem then is to select K such that J(R) is 
minimized. That is, 



K = argmin/(/?) 



(9) 



Let us define A = {A,- € 0,,j = 1,2, ...«} that is, instead of 
specifying an exact location for eigenvalues, the requirement 
is now to place the poles in a region defined as an union of 
disconnected closed regions on the complex plane. Assume 
each Q.i,i = 1,2,...,« are convex. Then the problem is to 
find the most optimum pole location and consequently the 
controller that yields minimum quadratic condition number 
of the closed loop modal matrix, i.e. 



A 



opt 



argminmax-jT^T 



(10) 



one can think of gradient based minimization techniques 
to find A pt, but in general the cost function, can be 
very arbitrary non-convex with respect to closed loop pole 
locations. 

This non-convex nature is illustrated in the subsequent 
example. 
Numerical Example 1: Consider a LTI system described as 



y(t) = CR[e D ^L T x(0) + / e D ^-^l} ' Br(x)dx\ (5) 
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Eigenvalues, together with eigenvectors plays an important 
role in shaping response. Judicial choice of eigenvector can 
prevent a mode from being excited by input. Using the 
control law u = ky, max(ra,r) eigenvalues may be assigned 
and min(m,r) entries of each corresponding eigenvector can 
be chosen precisely [3]. Assume for simplicity that the set 
of open and closed loop eigenvalues have no element in 
common), that is \Xj — A\ ^ for i = 1,2, ...n then 



Ri = (XiI-Ay l BKCRi = (XiI-Ay l BW, 



(6) 



where Wi = KCRj Eqn (6) can be written in matrix form 

\Wi W 2 ...W„] = KC[R X R 2 ...R„] (7) 



K 



W{CRY 



(8) 



The determination of controller matrix K therefore is 
straightforward. Note that the accuracy of computing K is 
actually depends on condition number of right eigenvector 
matrix. 




The desired closed loop pole locations are Ai = — 4.5,A2 = 
—6.5, A3 G [—5,0], and A4 = —3. The cost function 
J(R,X) = \\R\\ 2 \\R~ l \\ 2 where R is the closed loop modal 



matrix. To construct an approximate nature of cost function 
variation with A3, we took 500 samples of A3 variation. 
In each case we use the Rank-n Update eigenstructure 
assignment method [16] to find the minimum value of 
J(R*,X?, = Xit,i = 1,2,. ..500), and each case the method 
has converged giving best eigenstructure for that choice of 
A3. The figure clearly showing a non convex nature of the 
cost function. In general the nature of the cost function can 
be very arbitrary and gradient based minimization approach 
may lead to local minima. One way to avoid the problem 
is to start gradient based minimization with several initial 
condition, but in general these are difficult to estimate 
a priori. Another viable approach is to use Monte Carlo 
technique for finding probably approximate minima. In 
Monte Carlo procedure, instead of finding true J* one tries 
to find an estimate of J* with an accuracy e and confidence 
1 — 5. From classical one sided Chernoff bound one can 
find the minimum number of i.i.d samples to estimate the 
quantity J* with confidence 1 — 5 and accuracy e. More 
specifically if M > ^ln^ then Pr{\J* -J*\ < e} > 1 - 8. 
Now, we present a very simple Monte Carlo (non adaptive 
random search method) based procedure for low sensitivity 
eigenstructure assignment. 

Algorithm: Low Sensitivity Eigenstructure Assignment 
Generate M i.i.d self conjugate pole location set 5(A,-) 
according to a probability distribution (P x ) 

For i <— 1 to M do 

calculate min 77 and K for S'(Aj) using any method 

> One can use Direct method/ Rank N update/ Rank 

> one Update/ KNV method etc. 

Return The empirical optimum pole locations 
S{Xi) = arg min rj*(S(A,-)) 

i 

Below is an example of low sensitivity state feedback 
eigenstructure assignment using Monte carlo method. 
Numerical Example 2: Consider the previous LTI model 
defined in numerical example 1. The system is controllable. 
We designed a low sensitivity state feedback eigenstructure 
assignment for the following closed loop pole specifications. 
The closed loop pole location are specified as Ai € {( — 3.5 + 
/?cos0)+;(3+/?sin0),/f <= [O,],0e [0,-60]}, A 2 = Ai, 
A3 € [—1.6,-0.8] and A4 = —3. For an accuracy of 0.01 
and confidence level of 0.99 we need to generate at least 
23,025 samples to estimate the value of J*. We used a 
Rank-n-Update method to find the eigenvector corresponds to 
minimum quadratic condition number for a given eigenvalue. 
In this case, for every sampled eigenvalue in the specified 
region the Rank n Update algorithm converged. Finally we 
get Tl(R) = 3.6083. For simplicity we have not assume 
any structure of the eigenvector matrix. The corresponding 
controller 



and the eigenvalue locations are 



Ai = 


= -1.5783 + 2.4560/ 


A 2 = 


= -1.5783-2.4560* 


A 3 = 


= -0.9158 


A4 = 


= -3.0000 



K 



0.2188 
-5.2205 



0.0734 
-0.7562 



-0.7441 
-0.6763 



-0.5639 

3.3718 



IV. Robust Control Design using Eigenstructure 
Assignment 

Modelling errors are inevitably introduced when a model 
is obtained by approximating the system dynamics. These 
errors can be introduced in the form of structured / parametric 
or unstructured uncertainty. Due to these uncertainties, the 
eigenvalue of nominal system and the true system may differ, 
resulting in a difference in response characteristics. Low 
sensitivity eigenstructure assignment is an indirect method 
to tackle the effect of uncertainty but as shown in [17] that 
a small improvement in the sensitivity measure may result 
in a significant degradation in frequency domain measure 
of robustness. In contrast to the low sensitivity eigenstruc- 
ture assignment, robust control design by H^ optimization 
methods try to find a feedback controller that guarantees 
closed loop stability and performance robustness under all 
uncertainties. But the price to pay for this is, it often 
compromises with system bandwidth, as shown in [22], re- 
sulting sluggish response characteristics. Since time domain 
response characteristics governed by system eigenstructure, 
are equally important design concern, therefore it is useful to 
opt for a design method which satisfies both the time domain 
performance specifications provided by eigenstructure and 
frequency domain robust stability and performance specifi- 
cation. For MIMO systems, for a specific pole assignment 
problem, resulting controller is non unique. Therefore one 
can select a controller gain K which optimizes a Robust 
performance index. Early contribution on this area was by [6] 
and [18]. These works are concerned with optimization of a 
mixed sensitivity performance index as robustness measure 
under unstructured uncertainty. In the present work we tried 
to investigate the same for structured uncertainty. 

First consider the problem of robust stability using eigen- 
structure assignment under structured uncertainty. Earlier, 
Sobel et al [19] tried to solve this problem by placing the 
closed loop eigenvalues far away form the jco axis. 

Although, only eigenvalue assignment can guarantee ro- 
bust stability, but it alone cannot guarantee robust perfor- 
mance. Moreover the robust stability result is very conser- 
vative and not many design freedom exists. Therefore it is 
worth to look at the design freedom offered by simultaneous 
eigenvalue and eigenvector assignment. Let us formulate the 
problem as to ensure stability and performance robustness 
under all perturbation by designing suitable eigenstructure. 
Denote uncertainty affecting the LTI system described as D = 
{8 e F : 8 = bdiag(qiI u q 2 h,--q,JmAu—^b)} and B D = 
{8 eV>: \\q\\ p <p,o(Ai) <p,i = 1,2,. ..b}. We consider the 
family of uncertain systems obtained when the uncertainty 



8 varies over Bjj(p) and we say that a certain property, e.g. 
stability or performance is robustly satisfied if it satisfied 
for all members of the family. The problem becomes: For 
a plant with such structured/parametric uncertainty design a 
controller using eigenstructure assignment method so that the 
following robust stability and performance requirements are 
satisfied for all members in the family of uncertain plants. 



they have finite VC dimension and thereby poses UCEM 
property. However sample complexity bounds in this case 
are two large to be useful. Hence, we choose a more general 
class of algorithm to design a probabilistic eigenstructure 
assignment. 

Algorithm: Probabilistic Eigenstructure Assignment 



1) for all X p , Re(Xp) < i.e. perturbed plant will be Choose integers M > 



In i and N > MA In W- Gen 



robustly stable 

f 5 1 
2) sup a _ < e where S(s) and T(s) are the sensitivity 

and complementary sensitivity transfer function. 

During recent years it has been reported that several 
problem in robustness analysis and synthesis are either NP 
complete or NP hard. A survey of NP hard problems in con- 
trol theory can be found in the work of Blondel and Tsitkilis 
[7], Nemirovskii [20], Poljak and Rohn [21]. These results 
shows that various problems in robust control are practically 
unsolvable if the number of variables become sufficiently 
large. The problem of checking robust stability and per- 
formance under structured/ parametric uncertainty is proven 
as NP-Hard. However, parametric uncertainty can also be 
handled in H„ framework, but then the design will be too 
conservative, jj, Analysis/Synthesis approach although can 
handle parametric uncertainty, but computational complexity 
of ji calculation is again NP Hard. Also stability evaluation 
of interval matrices is itself a NP hard problem. To overcome 
the intractability issues and conservatives, recent trends are to 
consider the problem of stability and performance robustness 
in probabilistic framework [24]. Define each triple (A 7 B,C) 
where (A,B,C) are system matrices, belongs to Ax B xC. 
Denote A x B x C as S. Now suppose P is a probability mea- 
sure on the set S. Define e be a given accuracy parameter. For 
each matrix K eK define T(K) = {{A,B,C) : J t (A,B,C) < 
Yi,i = 1,2...«}. In probabilistic framework the question is: 
Does there exist a matrix K eK such that P(T(K)) > 1 - e? 
In other words, instead of looking for a matrix K e K which 
satisfies worst case performance, one compromises the search 
for a K e IK satisfies the performance index for most of the 
cases except possibly those belonging to set of volume no 
larger than e. 

Now we describe one important notion of probable ap- 
proximate minima as found in [14]. 

Definition: Suppose / : Y — ► R, that Py is a given prob- 
ability measure on Y and that e,a > are given numbers. 
A number /o e R is said to be a near minimum of /(.) to 
accuracy e and level a, or a probably approximate minimum 
of /(.) to accuracy e and level a if /o > /* — e and in 
addition P Y {y G Y : f(y) < f - e} < a. 

Based on this definition sample complexity bound has 
been established in [14], and two algorithms have been 
proposed. The first one is universal while the second 
one is applicable to situations where associated family of 
functions has the uniform convergence of empirical mean 
(UCEM)property. Now, a collection of sets can have UCEM 
property only when its VC dimension [14] is finite. For the 
class of function T(K). In [13] it has been established that 



2? m TT 



erate N i.i.d samples xi,X2, ■••*# from set S according to (P x ) 
and M i.i.d samples k\,ki,...k m from K according to Pf, 
For ;' <— 1 to M do 

M 

calculate /; = ^ I T(kj,Xi) 

Return The empirical optimum controller 

K(kj) = argmin/;- 

j 

Then with confidence 1 — 5, it can be said that /o is a 
probably approximate near minimum of /(.) to accuracy e 
and level a. 

Numerical Example 3: Consider a LTI plant model de- 
scribed as 
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The perturbation model considered as A = diag{8n, c>22, fe} 
where 8 n G [-1,0],§22 G [-1,0],5 33 G [-1,0]. The controller 
has the structure K = WR^ 1 where R is the right eigenvector 
matrix and W is the design vector matrix. Let's assume we 
have right eigenvector structure constrained as 



R 



x 1 
1 x 
1 x 



where a x signifies arbitrary values are possible. Therefore 
we define K = K{yi,y2,y$) where yi,y2,y$ are the diagonal 
elements of right eigenvector matrix. Define cost function as 



./ 



1 if(A(x),B,/sT,C)unstable 
— — otherwise 



1 + \Ul-GK)-l\ 



selecting e = 0.01, a = 0.01,5 = 0.05 results in M > 367 
and N > 22339. Desired pole locations are —1,-2,-5. After 
running the algorithm the controller matrix found as, 



K 



-3.0566 
3.1167 



-0.1954 
3.5937 



1.2314 
1.5775 



The right eigenvector matrix 



R 



-0.2336 0.5461 0.0927 

0.1701 -0.5118 0.0556 

-0.2766 0.8865 0.6306 



and corresponding ||5(s)||oo = 3.54. Note that this controller 
does not stabilizes the plant under all uncertainties. As a 
counter example, if the perturbed matrix is 

-0.9002 2 

A(A)= -1.8709 2 

-2 -2 -3.4576 



then the closed loop system has eigenvalue 
{0.033,-3.76,-5.5}. 

V. Conclusion 

We have shown how randomized algorithm can be used 
in Low sensitivity eigenstructure assignment for relaxed 
pole locations and robust control design using eigenstructure 
assignment method under structured perturbation. The al- 
gorithms are computationally straightforward to implement, 
although the estimate of sample complexity is conservative. 
Numerical examples are shown for state feedback eigenstruc- 
ture assignment although these procedure can be extended for 
output feedback. 
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